Numerical study of the noninertial systems: applicationto train coupler systems

Car coupler forces have a significant effect on the longitudinal train dynamics and stability. Because the coupler inertia is relatively small in comparison with the car inertia; the high stiffness associated with the coupler components can lead to high frequencies that adversely impact the computational efficiency of train models. The objective of this investigation is to study the effect of the coupler inertia on the train dynamics and on the computational efficiency as measured by the simulation time. To this end, two different models are developed for the car couplers; one model, called the inertial coupler model, includes the effect of the coupler inertia, while in the other model, called the noninertial model, the effect of the coupler inertia is neglected. Both inertial and noninertial coupler models used in this investigation are assumed to have the same coupler kinematic degrees of freedom that capture geometric nonlinearities and allow for the relative translation of the draft gears and end of car cushioning (EOC) devices as well as the relative rotation of the coupler shank. In both models, the coupler kinematic equations are expressed in terms of the car body and coupler coordinates. Both the inertial and noninertial models used in this study lead to a system of differential and algebraic equations that are solved simultaneously in order to determine the coordinates of the cars and couplers. In the case of the inertial model, the coupler kinematics is described using the absolute Cartesian coordinates, and the algebraic equations describe the kinematic constraints imposed on the motion of the system. In this case of the inertial model, the constraint equations are satisfied at the position, velocity, and acceleration levels. In the case of the noninertial model, the equations of motion are developed using the relative joint coordinates, thereby eliminating systematically the algebraic equations that represent the kinematic constraints. A quasistatic force analysis is used to determine a set of coupler nonlinear force algebraic equations for a given car configuration. These nonlinear force algebraic equations are solved iteratively to determine the coupler noninertial coordinates which enter into the formulation of the equations of motion of the train cars. The results obtained in this study showed that the neglect of the coupler inertia eliminates high frequency oscillations that can negatively impact the computational efficiency. The effect of these high frequencies that are attributed to the coupler inertia on the simulation time is examined using frequency and eigenvalue analyses. While the neglect of the coupler inertia leads, as demonstrated in this investigation, to a much more efficient model, the results obtained using the inertial and noninertial coupler models show good agreement, demonstrating that the coupler inertia can be neglected without having an adverse effect on the accuracy of the solutio

Natural Radioactivity in Some Building Materials Originating from a High Background Radiation Area

Twenty four samples of building materials, collected from utilized quarries dispersed randomly in a high natural background radiation area, were analyzed for 226Ra, 232Th and 40K by γ-spectrometry. This area lies in Hail province, Saudi Arabia. The collected samples were fragmented granites, granite gravels with clays, sands and crushed black rocks (mafic metavolcanic rocks). The results showed that the highest activity concentrations were found in the fragmented granite materials and ranged from 144-207, 671-1058 and 964-1440 Bq/kg with average values of 194, 912 and 1320 Bq/kg for 226Ra, 232Th and 40K, respectively. The lowest activity concentrations were found in the black rock materials which ranged from 19-39, 47-125 and 212-306 Bq/kg with average values of 24, 82 and 255 Bq/kg for 226Ra, 232Th and 40K, respectively. The radioactivity levels in the other materials lie somewhere in between. Granites and clays exceeded the proposed hazard indices for the usage as building materials and should be restricted, whereas the sands and the crushed black rocks complied with these indices and can be used without restrictions

Coupling dynamics of a geared multibody system supported by Elastohydrodynamic lubricated cylindrical joints

A comprehensive computational methodology to study the coupling dynamics of a geared multibody system supported by ElastoHydroDynamic (EHD) lubricated cylindrical joints is proposed throughout this work. The geared multibody system is described by using the Absolute-Coordinate-Based (ACB) method that combines the Natural Coordinate Formulation (NCF) describing rigid bodies and the Absolute Nodal Coordinate Formulation (ANCF) characterizing the flexible bodies. Based on the finite-short bearing approach, the EHD lubrication condition for the cylindrical joints supporting the geared system is considered here. The lubrication forces developed at the cylindrical joints are obtained by solving the Reynolds’ equation via the finite difference method. For the evaluation of the normal contact forces of gear pair along the Line Of Action (LOA), the time-varying mesh stiffness, mesh damping and Static Transmission Error (STE) are utilized. The time-varying mesh stiffness is calculated by using the Chaari’s methodology. The forces of sliding friction along the Off-Line-Of-Action (OLOA) are computed by using the Coulomb friction models with a time-varying coefficient of friction under the EHD lubrication condition of gear teeth. Finally, two numerical examples of application are presented to demonstrate and validate the proposed methodology.National Natural Science Foundations of China under Grant 11290151, 11221202 and 11002022, Beijing Higher Education Young Elite Teacher Project under Grant YETP1201

A study on wear evaluation of railway wheels based on multibody dynamics and wear computation

The wear evolution of railway wheels is a very important issue in railway engineering. In the past, the reprofiling intervals of railway vehicle steel wheels have been scheduled according to designers' experience. Today, more reliable and accurate tools in predicting wheel wear evolution and wheelset lifetime can be used in order to achieve economical and safety benefits. In this work, a computational tool that is able to predict the evolution of the wheel profiles for a given railway system, as a function of the distance run, is presented. The strategy adopted consists of using a commercial multibody software to study the railway dynamic problem and a purpose-built code for managing its pre- and post-processing data in order to compute the wear. The tool is applied here to realistic operation scenarios in order to assess the effect of some service conditions on the wheel wear progression

A comparison of Finite Elements for Nonlinear Beams: The absolute nodal coordinate and geometrically exact formulations

Two of the most popular finite element formulations for solving nonlinear beams are the absolute nodal coordinate and the geometrically exact approaches. Both can be applied to problems with very large deformations and strains, but they differ substantially at the continuous and the discrete levels. In addition, implementation and run-time computational costs also vary significantly. In the current work, we summarize the main features of the two formulations, highlighting their differences and similarities, and perform numerical benchmarks to assess their accuracy and robustness. The article concludes with recommendations for the choice of one formulation over the other

Modeling of the condyle elements within a biomechanical knee model

The development of a computational multibody knee model able to capture some of the fundamental properties of the human knee articulation is presented. This desideratum is reached by including the kinetics of the real knee articulation. The research question is whether an accurate modeling of the condyle contact in the knee will lead to reproduction of the complex combination of flexion/extension, abduction/adduction and tibial rotation ob-served in the real knee? The model is composed by two anatomic segments, the tibia and the femur, whose characteristics are functions of the geometric and anatomic properties of the real bones. The biomechanical model characterization is developed under the framework of multibody systems methodologies using Cartesian coordinates. The type of approach used in the proposed knee model is the joint surface contact conditions between ellipsoids, represent-ing the two femoral condyles, and points, representing the tibial plateau and the menisci. These elements are closely fitted to the actual knee geometry. This task is undertaken by con-sidering a parameter optimization process to replicate experimental data published in the lit-erature, namely that by Lafortune and his co-workers in 1992. Then, kinematic data in the form of flexion/extension patterns are imposed on the model corresponding to the stance phase of the human gait. From the results obtained, by performing several computational simulations, it can be observed that the knee model approximates the average secondary mo-tion patterns observed in the literature. Because the literature reports considerable inter-individual differences in the secondary motion patterns, the knee model presented here is also used to check whether it is possible to reproduce the observed differences with reasonable variations of bone shape parameters. This task is accomplished by a parameter study, in which the main variables that define the geometry of condyles are taken into account. It was observed that the data reveal a difference in secondary kinematics of the knee in flexion ver-sus extension. The likely explanation for this fact is the elastic component of the secondary motions created by the combination of joint forces and soft tissue deformations. The proposed knee model is, therefore, used to investigate whether this observed behavior can be explained by reasonable elastic deformations of the points representing the menisci in the model.Fundação para a Ciência e a Tecnologia (FCT) - PROPAFE – Design and Development of a Patello-Femoral Prosthesis (PTDC/EME-PME/67687/2006), DACHOR - Multibody Dynamics and Control of Hybrid Active Orthoses MIT-Pt/BSHHMS/0042/2008, BIOJOINTS - Development of advanced biological joint models for human locomotion biomechanics (PTDC/EME-PME/099764/2008)

Spatial rigid-multi-body systems with lubricated spherical clearance joints : modeling and simulation

The dynamic modeling and simulation of spatial rigid-multi-body systems with lubricated spherical joints is the main purpose of the present work. This issue is of paramount importance in the analysis and design of realistic multibody mechanical systems undergoing spatial motion. When the spherical clearance joint is modeled as dry contact; i.e., when there is no lubricant between the mechanical elements which constitute the joint, a body-to-body (typically metal-to-metal) contact takes place. The joint reaction forces in this case are evaluated through a Hertzian-based contact law. A hysteretic damping factor is included in the dry contact force model to account for the energy dissipation during the contact process. The presence of a fluid lubricant avoids the direct metal-to-metal contact. In this situation, the squeeze film action, due to the relative approaching motion between the mechanical joint elements, is considered utilizing the lubrication theory associated with the spherical bearings. In both cases, the intra-joint reaction forces are evaluated as functions of the geometrical, kinematical and physical characteristics of the spherical joint. These forces are then incorporated into a standard formulation of the system’s governing equations of motion as generalized external forces. A spatial four bar mechanism that includes a spherical clearance joint is considered here as example. The computational simulations are carried out with and without the fluid lubricant, and the results are compared with those obtained when the system is modeled with perfect joints only. From the general results it is observed that the system’s performance with lubricant effect presents fewer peaks in the kinematic and dynamic outputs, when compared with those from the dry contact joint model.Fundação para a Ciência e a Tecnologia (FCT

Comparison and Implementation of a Rigid and a Flexible Multibody Planetary Gearbox Model

We propose algorithms for developing (1) a rigid (constrained) and (2) a flexible planetary gearbox model. The two methods are compared against each other and advantages/disadvantages of each method are discussed. The rigid model (1) has gear tooth reaction forces expressed by Lagrange multipliers. The flexible approach (2) is being compared with the gear tooth forces from the rigid approach, first without damping and second the influence of damping is examined. Variable stiffness as a function of base circle arc length is implemented in the flexible approach such that it handles the realistic switch between one and two gear teeth in mesh. The final results are from modelling the planetary gearbox in a 500 kW wind turbine which we also described in Jørgensen et.al (2013)